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We study spin relaxation in dilute magnetic semicondnctors near a ferromagnetic transition, where 
spin fluctuations become strong. An enhancement in the scattering rate of itinerant carriers from 
the spin fluctuations of localized impurities leads to a change in the dominant spin relaxation mecha¬ 
nism from Dyakonov-Perel to spin flips in scattering. On the ferromagnetic side of the transition, we 
show that due to the presence of two magnetic components - the itinerant carriers and the magnetic 
impurities - with different gyromagnetic ratios, the relaxation rate of the total magnetization can 
be quite different from the relaxation rate of the spin. Following a disturbance of the equilibrium 
magnetization, the spin is initially redistributed between the two components to restore the equilib¬ 
rium magnetization. It is only on a longer time scale, controlled by the spin-orbit interaction, that 
the total spin itself relaxes to its equilibrium state. 


I. INTRODUCTION 

The draw of magnetic semiconductors is readily appar¬ 
ent from the name; they offer an opportunity to blend 
magnetic and semiconducting properties in a single ma¬ 
terial. In spintronics, this is especially relevant because of 
the additional means of interaction through the magnetic 
and electrical properties of carriers. Dilute magnetic 
semiconductors are unique in that they marry a low- 
magnetism ferromagnet with a semiconductor capable of 
maintaining long spin lifetimes. This makes possible the 
manipulation of carrier spins by ferromagnetic switching 
or, in the non-magnetic state, through more traditional 
means like optical excitation. While these materials are 
a reality, the ferromagnetic transitions are so far confined 
to well below room temperature.Nonetheless, the po¬ 
tential for these materials warrants their study, not to 
mention their suitability for studying spin systems near 
critical transitions, like a ferromagnetic transition. 

In this paper, we examine the impact of a ferromag¬ 
netic transition on carrier spin lifetimes. Strong spin fluc¬ 
tuations are expected for the constituents driving a fer¬ 
romagnetic phase change.® These fluctuations can man¬ 
ifest themselves through enhanced carrier scattering. In 
the Dyakonov-Perel spin relaxation mechanism,^® spin 
lifetimes are typically inversely related to scattering life¬ 
times. This is because the relaxation process occurs be¬ 
tween collisions, due to precession of the spin in the spin- 
orbit held. Then, in the Dyakonov-Perel mechanism, en¬ 
hanced scattering due to spin huctuations disrupts spin 
precession and can inhibit the rate of spin relaxation. 
Under these circumstances, we can hnd a crossover in 
the dominant spin relaxation mechanism to one based on 
spin hips in scattering. Spin hips can occur in scattering, 
for example, due to the mixing of spin-up and spin-down 
states from a spin-orbit interaction; this is the premise 
of the Elliott-Yafet spin relaxation mechanism. ^^’^2 
ferromagnetic system, an even greater rate of spin hips in 
interaction is made possible by transfer of spin between 


itinerant carriers and magnetic impurities. 

We model a ferromagnetic transition in a dilute mag¬ 
netic semiconductor with a degenerate gas of itinerant 
carriers (electrons or holes) interacting via spin exchange 
with dilute localized magnetic impurities. GaMnAs is 
the prototypical material to which this model has been 
applied - not without controversy over the hole transport 

mechanism. ^d3.14 QgGdN jg another example: here it is 

the s-f exchange interaction (as opposed to s-d interac¬ 
tion) that drives the ferromagnetic transition. ^®d5,i6 

this paper, we will keep the discussion fairly general, ap¬ 
plying parameters related to a GaMnAs system for the 
sake of calculations and analysis, but noting that the 
model is meant more as an example of critical phenom¬ 
ena than an exact description of GaMnAs. Other arti¬ 
cles have examined the critical behavior of GaMnAs in 
more detail.Our treatment is mean-field-like in 
that spin fluctuations and the interaction they mediate 
between the carriers are at the level of a random phase 
approximation (RPA). Thus, we ignore critical fluctua¬ 
tions which presumably lead to large corrections to the 
quasiparticle spectrum in the immediate vicinity of the 
ferromagnetic transition. 

Derivations of spin relaxation times resulting from spin 
dependent scattering are presented for two mechanisms: 
transfer of spin from carriers to magnetic impurities and 
Dyakonov-Perel. The results we obtain are applicable to 
dynamic, spin-polarized systems and focus primarily on 
the contribution of carrier-carrier interactions mediated 
by the spin fluctuations of magnetic impurities. We ap¬ 
ply common Fermi liquid theory techniques, beginning 
with the Kadanoff-Baym kinetic equation. By utiliz¬ 
ing the GW approximation for the carrier self-energy, 
we obtain relatively simple and easy to calculate ex¬ 
pressions for the spin relaxation rate, which are gener¬ 
ally applicable to a wide variety of effective interactions. 
The results derived here follow previous analytic^®’®®^^® 
and computationaP®^^® studies, but extend them signif¬ 
icantly by including dynamics and spin exchange. 

This paper is organized as follows: in Section II, we 
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present an effective interaction that can be used to gen¬ 
erate a ferromagnetic transition and which will be used 
to calculate carrier scattering rates; in Section III, we 
discuss the ferromagnetic phase change and the behav¬ 
ior of the scattering amplitude near the transition; in 
Section IV we derive analytic expressions for the rate of 
spin relaxation due to spin-flips in scattering and the 
Dyakonov-Perel mechanism; in Section V, we discuss 
some results of spin relaxation across a ferromagnetic 
transition; Section VI contains our concluding remarks. 


II. EFFECTIVE SPIN-SPIN INTERACTION 

We first introduce the model dilute magnetic semicon¬ 
ductor on which we base calculations. Our model resem¬ 
bles the Zener model of ferromagnetism, with a ferro¬ 
magnetic transition being driven by spin-exchange inter¬ 
actions between itinerant carriers and localized magnetic 
impurities.The Hamiltonian includes contributions 
from carriers Ha, magnetic impurities "Hs, and interac¬ 
tions between carriers and magnetic impurities "Ho-s: 


^ ^ / 2 ^* Bcr(rz)'i , 

i ^ 

(la) 

Bs = g-sgs ^ Sj ■ Bs(Rj), 

(lb) 

BaS = - R-j)> 

(Ic) 




where and (t^ are, respectively, the positions and spins 
of the mobile carriers with effective mass m* , Rj and Sj 
are the positions and spins of the fixed impurities, /is is 
the Bohr magneton, and gs are the g-factors which 
determine the magnetic moments of the carriers and the 
magnetic impurities as gsga and gsgs, respectively, and 
Bg. and B 5 are magnetic fields which are supposed to 
act only on carriers or magnetic impurities, respectively. 
The contact-type interaction between carriers and mag¬ 
netic impurities is represented by J(5(ri —R^), where J is 
the interaction strength averaged over the unit cell. The 
impurities are taken to be sufficiently dilute that they do 
not interact with each other directly. 

The alignment of spins with an external magnetic field 
B is linearized as gsgi^) = x(B), where x is the mag¬ 
netic susceptibility. By combining the interaction J with 
the magnetic fields, we can derive effective mean fields 
B®® and Bg® acting on carriers and impurities, respec¬ 
tively: 

HBgaBfir) = gsgaB^ir) + JS(r), (2) 

Mb5sB5®(R) = gsgs^si^) + Jcr(R), (3) 


relatively small strength compared to J. The magnetic 
fields in Eqs. (2) and (3) are specific to each species, so 
it is cleanest to absorb gLsg into these fields and write 

(^)=xW(Bf), (4) 

(S)=xi(Bf), (5) 


where x^<y nnd x^ss non-interacting (with respect to 
J) spin-spin susceptibilities for carriers and magnetic im¬ 
purities, respectively. They have dimensions of inverse 
energy-volume. To be explicit, these fields relate to the 
physical magnetic field by Bs = gsfis^ and the suscep¬ 
tibilities are related to textbook magnetic susceptibilities 
by xss = xK.gg'sf- 

The equilibrium spin polarization of the impurities 
(S'z) follows the Brillouin function Bj{x),^° 


(5z) = -SsnsBss 


Ss{Bs,z) 

ksT 


( 6 ) 


where S 5 is the magnitude of the impurity spin, ns is the 
density of magnetic impurities, and T is the temperature. 
Then, the longitudinal and transverse components of the 
spin susceptibility of the impurities are calculated from 

xP^s.=d{Sz^)/d{Bs,z), (7) 

xT^s^=2{S.)/{Bs,z). ( 8 ) 

The non-interacting spin response of carriers coincides 
with Lindhard response: 


X, 


(0) 


(9,w) = 




koL 


k.a 


kio + Sko 


fO 

J k+qo! 


^k+qa 


xtV(9>^)=4S 


2 


fO 


fct 


huJ + Ek'f 


fO 

J k+q^ 




( 9 ) 

( 10 ) 


where So- is the magnitude of the carrier spin, is the 
equilibrium Fermi-Dirac distribution, and x^-cr+ig^'-^) = 
[xi+cT_(g,-w)]*. 

We can also shift the interaction J into the spin suscep¬ 
tibilities to generate effective response functions. This is 
accomplished by solving the coupled linear Eqs. (4) and 
(5) for (cr) and (S) and identifying the effective response 
function as the proportionality between average spin and 
bare field: 


■(5z)' 


{Bs,z) 

_(crz)_ 


{B.,z)_ 


and 


( 11 ) 


where cr(r) = <5(r-ri)cri and S(r) = Y,j <5(B--Bj)Sj 

r(^+)i 


'{Bs+y 

are the densities of carrier spin evaluated at r and im¬ 

(5_) 

1 

{Bs-) 

purity spin evaluated at R, respectively. We are not in¬ 

(c^+) 

“ 2 = 

{B.+) 

cluding spin-orbit in the effective carrier field due to its 



l{B.-)\ 
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FIG. 1. The effective spin-spin interaction between carriers. 
While carriers do not interact directly with each other via 
spin, the spin interaction J can be mediated by impurities 
through the effective spin-spin susceptibility xss- Here, a and 
a' are the incoming and outgoing carrier spins, respectively. 


where s± = Sx^isy and Bs± = Bs^x^iBg^y The effective 
susceptibilities are found to be 


and 

Xt = 


XS,S, XS.cr, 


Vxs°s, -J 

Xa^Sz X(7z<yz__ 


- 1 

Sb 

X 

I 

_I 


(13) 
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Rather than give explicit expressions for every effective 
response function in Eqs. (13) and (14), we write only 
those which are used in this paper. Namely, the effective 
impurity spin-spin response functions are 


Xs^sAq,^) = 


Xs+s- ( 9 ,w) = 


(0) 


1 - J‘^xs!s,Xa}aAq,^) 


(0) 

Xs+s. 


1 - {J/2AxT+S.Xala.{q,uj) 


(15) 

(16) 


where xs_s+(<?,w) = [xs+s_ (g,-w)]*. 

The interaction of interest to us is the effective carrier- 
carrier interaction, shown schematically in Fig. 1. By 
using the effective spin-spin susceptibility for impurities 
in Eqs. (15)-(16), we capture the effect of collective impu¬ 
rity spin fluctuations on carrier-carrier scattering. Notice 
that although the non-interacting impurity spin response 
is frequency independent, the effective spin response is 
dynamic due to the inclusion of carrier dynamics. The 
longitudinal (||) and transverse (T) components of this 
interaction are 


W^||q(w) = T^xs,s,,(g,w), (17) 

W'_Lq(a;) = (J/2)2 [xs+s_(g,w)-bxs_s+(g,w)] . (18) 


It can be verified that these effective interactions are 
equivalent to an infinite sum of all possible bare impurity- 
carrier interactions mediated by non-interacting suscep¬ 
tibilities, e.g. 


IE = ^ 

n—1 


( 0 )' 

Xss 



n—1 


III. FERROMAGNETIC TRANSITION AND 
PHASE 


The interaction J between carriers and impurities al¬ 
lows for the possibility of non-zero average spins (cr) and 
(S) in the absence of external fields. This is confirmed 
by the Curie-Weiss-type susceptibility to which the ef¬ 
fective response functions in Eqs. (15) and (16) reduce in 
the static, long wavelength limit: x (T — Tc)~^. The 
paramagnetic-ferromagnetic phase change occurs when 

1 - J^^Xs^s. (^c) lim xtS. (9> 0) = 0 

^ J"nsS2iV(0), (19) 

6kB 

where iV(0) is the carrier density of states at the Fermi 
level. 

We can self-consistently solve for the polarizations of, 
and average fields acting on magnetic impurities and car¬ 
riers below Tc- The polarization of the magnetic impuri¬ 
ties is given by the Brillouin function Bj{x): 


Ps 


Bss 


Ss\{BsA\ ' 

keT 


( 20 ) 


To a good approximation, the field acting on carriers is 
equivalent to the first order interaction with magnetic 
impurities: 


\{B„A\ = JSsnsPs- ( 21 ) 

The field acting on carriers also defines the spin splitting 
energy — £fi about the Fermi energy ep: 




£f 


(1 + P^)2/3_(1_P^)2/3 


( 22 ) 


Finally, we ensure self-consistency in the generated fields 
through expectation of a Goldstone mode on the ferro¬ 
magnetic side of the transition. Indeed, in this broken 
symmetry state the polarized spins can all rotate together 
without any energy cost. This mode is identified by set¬ 
ting g = a; = 0 in the effective transverse susceptibility 
and stipulating a divergence: 


1 - ( J/A'^Xsls.i'^) (g,o) = 


1-(J/2)2 


2{SA 2{aA 
{Bs,z) {Ba,z) 


0. (23) 
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T/T, T/T, 


FIG. 2. (Color online) Mean polarizations of, and effective 
magnetic fields acting on carriers and magnetic impurities 
have been calculated from simultaneous solution of Eqs. (20)- 
(23). GaMuAs parameters have been used to define the car¬ 
riers as holes with effective mass m* = O.Sme and density 
rio- = 5 X 10^° cm“®. The Mn acceptors have spin Ss = 5/2 
and density ns = 0.05/fIo, where flo = 0.45x 10”^^ cm® is the 
unit cell volume in GaAs. a) The carrier and impurity spins 
closely follow the approximation PsiT) = Ts(0)'\/l — {T/TcY 
from Ref. 31. The carrier polarization does not saturate, i.e. 
Po-(O) 7 ^ 1, due to a significantly smaller magnetic susceptibil¬ 
ity than for the impurities, b) As the effective fields are self¬ 
generated, their behavior parallels the ■y/r^^/T/T/)^ shape 
of the polarizations. The field acting on magnetic impurities 
is largely due to the second order RKKY-like interaction me¬ 
diated by carriers, and scales with for varying interaction 
strengths. 


The resulting polarizations and fields for several values 
of J are plotted in Fig. 2. 

It is instructive to examine the strength of the singu¬ 
larity in the effective interactions when Tc is approached. 
Approximating the non-interacting carrier susceptibility 
in the paramagnetic phase by its static, small-g form 
— ~S^A^(0) (l — /12ky), we can write the 

denominator of the scattering amplitude in Eq. (17) as 


We focus on two mechanisms of carrier spin relax¬ 
ation: transfer of spin to magnetic impurities (MI) and 
Dyakonov-Perel (DP). Under a wide variety of tempera¬ 
tures, scattering sources, and degeneracy, the DP mech¬ 
anism dominates in III-V semiconductors.®^’^® Thus, we 
find it relevant to include derivations for this spin relax¬ 
ation mechanism here. As will become evident, however, 
the inverse relationship between scattering and spin re¬ 
laxation in the DP mechanism hinders its effectiveness 
when the interaction that is responsible for carrier-carrier 
scattering becomes strong and long-ranged, e.g. near the 
ferromagnetic transition. This leads to a change in dom¬ 
inant carrier spin relaxation mechanism near the ferro¬ 
magnetic transition. Both the MI and the Elliott-Yafet 
(EY) mechanisms rely on spin flips in scattering; thus, 
they scale proportionally to the scattering rate and be¬ 
come more important when the carrier-carrier interac¬ 
tion becomes strong. This point is in agreement with 
Morandi, et who find that spin-flip-based spin re¬ 

laxation mechanisms (including a small Elliott-Yafet con¬ 
tribution) are important when spin-dependent scattering 
is dominant. However, the timescales of these two mech¬ 
anisms are very different. Since EY depends on the weak 
spin-orbit field to allow spin flips in interaction, it is gen¬ 
erally negligible when MI processes are occurring. This 
leads us to identify MI as the dominant spin relaxation 
mechanism near the ferromagnetic transition, and so we 
neglect the EY mechanism in the following derivations. 
Nevertheless, even when MI dominates, there is a role 
left for EY to further relax the total magnetization of 
impurities plus carriers after the total spin angular mo¬ 
mentum has reached equilibrium. This will be qualita¬ 
tively discussed in Section V. Notice that the enhanced 
carrier-carrier collision rate near the transition is not in 
opposition with, but rather the consequence of the crit¬ 
ical slowing down of the collective spin dynamics of the 
carrier-impurity system. 


A. Spin flips in scattering 


1 - ( 24 ) 

Then asT ^ T^, we have the same q~^ divergence in the 
interaction strength as a bare Coulomb interaction. 


IV. CARRIER SPIN RELAXATION 


We calculate the rate of spin relaxation due to spin 
flips in scattering using Keldysh formalism.®®’®^ The 
time rate of change of the average spin density (cTk) = 
Tr[/3k(T2] is 


dt 


ricrScr Tr 



(25) 


Typically, spin relaxation is understood to be the evo¬ 
lution of out-of-equilibrium polarized spins into an un¬ 
polarized equilibrium state. In a ferromagnetic system 
where there is an equilibrium spin polarization, spin re¬ 
laxation refers to the evolution of non-equilibrium spins 
into the polarized equilibrium state. Thus, rather than 
picturing the process as decay of polarization, we under¬ 
stand it as restoration of equilibrium polarization. 


where pk is a 2 x 2 density matrix in spin space for car¬ 
riers, (Tz is the z-Pauli matrix, and z has been chosen for 
the direction of spin polarization. Taking advantage of 
the relaxation time approximation for d{(7k)/dt and av¬ 
eraging over all wavevectors k, the spin relaxation time 
Ts is calculated from 


1 ^ Ek Tr [dtPk^z] 
Ts Ek Tr [pj^dz] 


( 26 ) 
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where is the first order non-equilibrium correction in 
the linearized density matrix. In the absence of spin- 
orbit and other external fields, the density matrix evolves 
according to 


_ f 

dt 


(27) 


where Ik is the 2 x 2 collision integral in spin space for 
carriers. The collision integral is expressed in terms of 
carrier self-energies Sk(w) and Green’s functions Gklw) 
as follows:^^’^^ 


dk 



-{s>(cc),G<(a;)}], (28) 


where {• ,•} represents an anti-commutation. In the ab¬ 
sence of spin-orbit interactions, small deviations from the 
equilibrium spin polarization only show up on the diag¬ 
onal elements of the density matrix. Then, we know a 
priori that Gk(w) and Sk(a;) are diagonal and can write 
the collision integral per spin-a as 

4. = / ^ - SL(^)G^Ja.)] . (29) 

We evaluate the self-energy in the GW 

approximation, 

(30) 

In allowing for arbitrary spin polarization, we sepa¬ 
rate the GW spin-space operator into longitudinal and 
transverse components according to Fig. 3: GW = 
W||iTzG(Tz -I- Wj_(i 7 a;G(Ta; -|- ayG&y). This allows us to 
rewrite Eq. (30) as 




+2Wfq(F!)G5k_q Ja; - fl)j , (31) 


where the newly defined longitudinal and transverse 
Green’s functions are Gyk = dzGkdz and G_Lk = 
(^/^) [^^cGkdz; -f (JyGkdy]. 

The GKB ansatz^^ - an exact relation in non¬ 
interacting systems - is utilized to relate the lesser and 
greater Green’s functions to the density matrices = pk 
and p^ = 1 — pk, where 1 is the 2 x 2 identity matrix: 

= TGU^)pt ± (32) 

Notice that we have neglected retardation effects in the 
density matrix by making it frequency independent. This 
is commonly referred to as the Markovian approximation 
and is justified when the time evolution of the distribu¬ 
tion function is slow on the time scale of microscopic 



FIG. 3. The GW self-energy includes separate contributions 
from longitudinal Wy and transverse Wj_ interactions. Here, 
a, /3, 7 , and 5 are spin states on which the interactions aiW&i 
operate. 


collisions. We let all non-equilibrium properties reside 
in pkj so G]f'{uS) are taken to be equilibrium Green’s 
functions. Furthermore, given that our aim is to derive 
a Boltzmann-like collision integral, we retain only the 
quasiparticle singularities in G^’“(a;), as those are the 
parts which conserve energy in collisions.^® Under these 
circumstances, the lesser and greater Green’s functions 
become 

GkQ(^) = ±*27r(5(;ia; - efca)Pka- (33) 

Similarly, we can use the fluctuation dissipation theorem 
to relate the lesser and greater interactions W^ (w) to the 
imaginary part of the retarded interaction 3mWq(u;):®® 

W^(w) = ±i2nB(±a;)9mWq(a;), (34) 

where nsiuj) = — 1)“^ is the Bose-Einstein distri¬ 

bution and P = \kBT)~^. 

The result of inserting Eqs. (31), (33), and (34) into 
the collision integral, Eq. (29), is 

/ka — 2 ^ ( / doj viTuIUljq (u^)S^k—qa IlCu) 

q 

[^i?(^)(f Pka)Pk—qa T ^i?( ^)7’ka(l Pk—qa)] 

-b 2'^mW±q{uj)S{eka - Ek-qa - ^) 

[^i?(a^)(l Pka)Pk—qa “b '^b{, a^)Pka(l Pk—qa)] 

(35) 

For a slightly out of equilibrium spin polarization, the 
distribution pka can be linearized as follows: 

Pka = /L + a/deafkM - /fea). (36) 

where £« is a small spin perturbation energy, a = -1-1 for 
spin-up, and a = — 1 for spin-down. Notice that the shifts 
Ea and Sa are related by the condition Eafkai^ ~ fka) = 
Ea/fea(l ~ /fea)I equivalently after summing over k- 
space, EaNaiO) = EaNaiO), meaning that the departure 
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a) J — 45 meV nm^ b) J — 60 meV nin^ c) J — 75 meV nrn^ 



dashed line in Fig. 4 for three different values of the ex¬ 
change coupling J, above and below the ferromagnetic 
transition temperature. The sharp drop of in the im¬ 
mediate vicinity of the ferromagnetic transition reflects 
the enhanced spin-flip rate from critical fluctuation in 
this region. 


B. Dyakonov-Perel mechanism 


FIG. 4. The spin-flip based spin relaxation time (solid lines) 
resulting from the spin dependent interaction is plotted vs. 
temperature for several values of J. For comparison, the 
Dyakonov-Perel spin relaxation time is shown by the dashed 
lines. We used GaAs parameters for holes in calculations: 
m* = 0.5me, cta.o. = 0.34, and Eg = 1.5 eV. The vertical 
lines mark the critical temperature in each plot. The short 
spin relaxation time near the ferromagnetic transition corre¬ 
lates with short spin lifetimes in the spin-flip mechanisms and 
long spin lifetimes in the Dyakonov-Perel mechanism. Away 
from the ferromagnetic transition, we see Dyakonov-Perel spin 
relaxation becoming dominant as is often found in III-V semi¬ 
conductors. 


from equilibrium involves only a change in spin density, 
not a change in carrier density. From a physical stand¬ 
point, it is clear that longitudinal interactions will pre¬ 
serve spin polarization, and so we should see that part 
of the collision integral vanish upon insertion of Eq. (36) 
into Eq. (35). Indeed, this is the case and we only need 
to evaluate 






fO _ £0 
J koL J k—qa: 

4 sinh^ {phoj /2) 


(37) 


which is entirely controlled by the transverse part of the 
interaction. 

Returning to the spin relaxation rate defined in 
Eq. (26), we can insert expressions for pk and /k. The 
rate of spin relaxation from transverse spin scattering 
events is 


1 

Ts 


1 

TT \N^{0) 


1 

mo) 


X 


f dui'^mW±q{uj) 

q,Q 


SwXga[l,g,tu] 
4sinh^(/3Iia;/2) ’ 


(38) 


where we have introduced the Lindhard-like imaginary 
response function 


TT d(ek-(-q/3 - Ska “ [/k-bq/3 “ fL] ^' (39) 

k 

The variable ^ is 1 in Eq. (38) but will take on more 
interesting forms in subsequent sections. The spin re¬ 
laxation time calculated from Eq. (38) is plotted as a 


In the DP mechanism, spins relax due to precession 
around axes defined by a magnetic field which varies in 
k-space. In III-V semiconductors, this field is typically 
Dresselhaus and/or Rashba spin-orbit. The effect of scat¬ 
tering on this process is usually to inhibit spin relaxation 
as spins are unable to make full precessions before be¬ 
ing relocated in momentum space. This sets a timescale, 
where the rate of scattering must be faster than the rate 
of spin precession for spin relaxation to be inhibited, oth¬ 
erwise spins enter a regime of free precessions and the rate 
of relaxation is directly proportional to the rate of pre¬ 
cession. In the limit of fast scattering, spins are unable to 
make appreciable precessions around the spin-orbit field’s 
axes and spin relaxation is entirely suppressed. This is 
the scenario we encounter. 

The derivation of DP spin relaxation is reviewed in 
Refs. 25, 39, and 40. When the Hamiltonian includes a 
contribution from a spin-orbit field Jlk of the form 

ni = Jok • <T, (40) 

where it is the Pauli matrix vector, the rate of spin re¬ 
laxation 1/ts for carriers due to the DP mechanism is^® 

1 2 Ek (^k) Tr [pgdz] 

3 EkTr[pOd,] ’ ^ ^ 


where r/ is an effective scattering time that weights col¬ 
lisions according to their ability to randomize the spin 
precession axis. The angle brackets (•) in Eq. (41) rep¬ 
resent an average over all directions of k, and in a bulk 
III-V semiconductor this quantity evaluates to 


i^l) 


16 o^L.gfc 
35 n?Eg ’ 


(42) 


where a^.o. is the spin-orbit coupling constant and Eg is 
the band gap energy. 

We emphasize here that spin-flip processes are not con¬ 
sidered in Eq. (41). This spin relaxation rate is entirely 
due to spin-orbit coupling. The challenge in evaluating 
Eq. (41) comes from the effective scattering rate I/t^. 
An oft-made approximation is to substitute the rate of 
momentum relaxation in its place. The validity of this 
approximation in the context of spin-scattering will be 
discussed in section V. In this section we present our cal¬ 
culation of 1 /t/. 
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Written in terms of density (/”) and spin (/'*) distri¬ 
butions, the quasi-equilibrium density matrix /5°. is 

Pk = /fe"i + (/fc* + (43) 

where 

fk = (fit - /fet)/ 2 > 

and was defined in Eq. (36). In writing the rate 
of spin relaxation in Eq. (41), there is an implicit as¬ 
sumption that a small spin polarization has been pre¬ 
pared which will relax. This is the reason for the non¬ 
equilibrium contribution in Eq. (43); it is stable on 
the short time scale of scattering, but variable on the 
longer time scale of spin relaxation. Then, while it is nec¬ 
essary to include this term in Eq. (41), it can be dropped 
when evaluating the collision integral, where we can in¬ 
stead use = /o„. 

The non-equilibrium part of the density matrix pik 
arises from interaction with the spin-orbit field:^® 

= (44) 

where [•, •] represents a commutation and & = 

(Ty, dz} is the Pauli spin vector. We use a relaxation 
time approximation to relate pj^ and to the collision 
integral 1^ as follows: 

-4=4. (45) 

One should exercise caution when evaluating this colli¬ 
sion integral, though. While Ik in Eq. (28) is still gen¬ 
erally applicable, the GKB ansatz for is not 

simply dependent on the difference ti — t 2 when a ^ P, 
and so does not lend itself easily to Fourier transform. 
Instead, we have (see Appendix A): 

GtpitiM) = (46) 

Thus, we opt to start in the time-domain with collision 
integral^"^ 


4(0= r dt' [i:<(t,t')G>{t',t) + G>[t,t')T.<{t',t) 

J —<50 

-E>(t,t')G<(t',0 - G<(t,t')S^(i',0] ■ (47) 

The self-energy in the GW approximation is 


sf(G,t 2 ) = [^lfq(G, 0 )G|_^(ti,t 2 ) 

q 


( 48 ) 


After inserting Eqs. (46) and (48) into Eq. (47), the time 
of observation can be set to t = 0 and integrals over 
t' performed. The fluctuation-dissipation theorem from 
Eq. (34) is again utilized for W^(w) and the collision 
integral reduces to a strictly off-diagonal matrix: 

/ oo 

dw |S5mW||q(w)(5(efcT, - ek-q 7 - 4;) 

([nf;( C^) -\- 4—q^] Pkaa [^i3(^) 4“ 47 ] Pk—qaa) 

- 29mlT_Lq(w)(5(efc7 - ek-q 7 - 4;) 


[^s( ^) 4” 4—q7] Pkaa^* 


(49) 


This is an appropriate time to re-introduce the relaxation 
time approximation from Eq. (45) and substitute the 
non-equilibrium density matrix elements from Eq. (44). 
Taking advantage of ea/L(l - 44) = ea4°a(l “ /4): 
the kinetic equation simplifies to 

i^£X(l-/fe7)^k+ = 

7 

/ duj< ^5mT^||q(a;)(5(£fc7 - ^k-q7 - 

„ ^ J-oo I 


q,7 


-pO _ pO 

Jk-v J] 


■'^7 ■'k-q7 + 4-qf^k-q+) 

^4sinh4/34;/2) ^ k q k q+; 

-I- 2 QmW±q{uj)S{skj - ek-q7 - 4;) 

fO _ fO 
J kj J k—q7 


Pe^ 


4 sinn (/34 j/2) 


4f2k+ 


(50) 


where r2k+ = ^k,y 4- fUk.x- 

This integral equation can be solved exactly using the 
methods of Sykes and Brooker,^^ but given the slowly 
varying nature of 4 around the Fermi level in the de¬ 
generate regime, it is a reasonable approximation to treat 
it as a constant evaluated at ep and extract it from the 
collision integral. Then, we only need to specify the spin- 
orbit field and solve for the effective scattering rate. 

The Dresselhaus field^^ components, written in terms 
of spherical harmonics are 


f^k,a; = 1/71/21 {Y3^i{d,(p) - Y3^_i{d,(p)) 

-f 47/735 (>3,3(79, ( fi ) - 4 - 3 ('d, ( fi )) 
f^k,y = %k^ (¥ 3 , 1 ( 79 , yj) -|-y3,_i(i9, 4) 

-71/7/735 (>3,3(79, ip) Y Ys^^sid, (/?)) 


(51a) 


(51b) 


Ilk,. = lloA:V84l05 (>^ 3.2 (t?, 4 + >3.-2(^^, 4), (51c) 


where Oq = otcf? j yTrofi//. These forms allow us to take 
advantage of the addition theorem for spherical harmon¬ 
ics, resulting in the following useful identity: 

J ^14-k'f^k',i = ^^k., J ^4-k'43(cOS79'), 

(52) 















where dVl' = d{cos'd')dip' and j can be any of x, y, or z. 
Pushing onto both sides of Eq. (50), we can 

extract the response function ‘^mxap j 9 ) '^] defined 


in Eq. (39), where 


^(dp) 


k 

3 

Pn 

rk.(k-bq)l 

Jk-f q|_ 

-^3 

. fc|k + q| _ 


( 53 ) 


Notice only the magnitude |q| is necessary to evaluate 
Q'wXa /3 ) 9 , w] since k is eventually integrated over all 
directions and can be measured relative to some arbitrary 
direction of q. The final result for the effective scattering 
rate due to a spin dependent interaction in the degenerate 
regime is 


1 _ ^ 

rf TT 

q.7 ' 


dw 


STOVE|!q(a;) q,u;] 2^771 V^Lg(a;) g,a;] 


NJO) 


4sinh^(/3/iw/2) 


N^(0) 4sinh^(/3?ia;/2)_ 


(54) 


The modified response function QmXap ) 9) can be evaluated analytically. It is rather cumbersome to write in 
full, so we provide the following form before the final integration is performed: 




TNpjO) 

49/3 




7riVc(0) 

49a 


0[1 




x-pp 

3/2 

Pq 

2 x-pp- qj 

X 

3 

2 ^x'^ - xpp 

X 

•ill 

Pi 

2 x + Pa- 9a 

X + Pa 

2y/a;2 -1- xpa 


(55) 


where 


qj = q/kp^, 

huj + {P — a)A/2 ^ q 

hqvpj 2kpj’ 

hio + (p — a) A/2 

di = -;-, 

epj 

and A is the same spin splitting defined in Eq. (22). Note that although a and P are not explicitly indexed in v^± 
and rj-y, their ordering in /3 — a is based on direct substitution of these quantities into Eq. (55). To be explicit, 
when evaluating 'Amx/3a, the order changes to a — p. For comparison, the effective scattering rate obtained when a 
spin-independent interaction Wq is used is 


J_ _ ^ f°° ^ ’AmWqju;) 'AmXj-y [1,9, ^] - ,q,uj] /spin-indep. \ 

^ ^ 7 ( 0 ) 4sinh^(/3Sa;/2) Vinteractiony 


(56) 


In Fig. 5 we plot the calculated values of t/r/^ vs. temperature for different values of exchange coupling and carrier 
concentration. We also plot the ordinary momentum relaxation rate which is related to the Drude resistivity Pap 
through = '^i^npe^pap/m*. A standard calculation yields 


1 


'^kp'a 


P 

7riV„(0) 



3mIT||q(w) 


tuxj + Eq SmXaa[l,g,a;] 

£Fa 4sinh^(/3Iia;/2) 


2’AmWpq{(jj) 


huj + Eq-t- a A 9mXaa[l, 9, w] 
spa Asi\ih^{phjj/2) 


We note that this formula generalizes a previous ex¬ 
pression that was derived for the contribution of spin fluc¬ 
tuations to the momentum relaxation rate (see Eq. (18) 
of Ref. 19). The essential difference is that our formula 
takes into account inelasticity of the scattering, whereas 


Ref. 19 treated the spin fluctuations in a quasi-static 
limit, which misses the critical enhancement of the scat¬ 
tering rate at the ferromagnetic transition. The quasi¬ 
static formula is obtained from the present formula by 
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a) J-dependence b) ricr-dependence 



FIG. 5. (Color online) The effective scattering rate used in 
Dyakonov-Perel spin relaxation calculations (solid lines) is 
compared to the momentum relaxation rate (dashed lines) for 
varying a) exchange interaction strengths and b) carrier densi¬ 
ties. The solid vertical lines mark critical temperatures. The 
effective scattering rates diverge at the ferromagnetic transi¬ 
tion, whereas the rate of momentum relaxation does not. This 
ultimately comes from zero-momentum transfer processes in 
spin-spin scattering events that affect Dyakonov-Perel spin 
relaxation, but not momentum relaxation. 


the replacements 


3mFF||(_L)q(w) -/3|;7||(_L)(q)|^a;(5(w), (57) 

where |hf||(_L)(q)p is the squared matrix element of the 
spin-dependent part of the static electron-impurity in¬ 
teraction (longitudinal or transverse) between the initial 
and the final state of the carrier. Note that one can¬ 
not set U!S(uj) ~ 0 in Eq. (57) because that distribution 
will be multiplied, in Eq. (57), by quantities that tend to 
infinity for w —>■ 0. 


V. DISCUSSION 

Our expression (54) for the effective scattering rate in 
the DP mechanism has some noteworthy features that 
make it qualitatively different from the previous known 
result in Eq. (56), which holds for spin-independent in¬ 
teractions. In particular, notice that there is a sign dif¬ 
ference in front of 3771 X 77 9 j'j-']- The significance 

of this sign is easier to understand if we look back to 
Eq. (50). In the longitudinal term, after extracting as 
a constant, we have a factor (UkT. -|-Dk_q_|.). This should 
be compared to (Dk+ — Dk-q+) which appears when a 
spin-independent interaction is used. The minus sign in 
the second case bears some similarity to the (1 — cosO) 
factor that appears in the standard expression for mo¬ 
mentum relaxation. A similar factor can be extracted 
from (Dk-(- — f^k-q+) and has the form [1 — P 3 (cos 0 )], 
where Paix) is the n = 3 Legendre polynomial.Eor 
the case of momentum relaxation, when there is no rela¬ 
tive angle 9 between incoming and outgoing momenta in 


a scattering event, there is no contribution to momentum 
relaxation. The analogy for the effective scattering rate 
is that zero-momentum transfers do not relocate particles 
in a k-dependent magnetic field, so spins can continue to 
process around the same magnetic field axis and those 
scattering events do not reduce the rate of spin relax¬ 
ation. This analogy is not appropriate for spin-dependent 
scattering, where processes with zero momentum trans¬ 
fer may still change spin orientation. In fact, what we 
find is a sort of generalized DP spin relaxation, where 
spin-dependent and spin-independent scattering should 
both be taken into account for the overall disruption of 
spin precessions in a spin-orbit field. 

The effective scattering rate from Eq. (54) is compared 
to the momentum relaxation rate from Eq. (57) in Eig. 5 
for varying interaction strengths and carrier densities. 
Unlike previous reports, the effective scattering rate is 
found to be generally larger than the momentum relax¬ 
ation rate.^^ This is consistent with the idea that more 
scattering events disrupt spin relaxation than affect mo¬ 
mentum relaxation. Especially near the critical temper¬ 
ature Tc, notice that momentum relaxation has a finite 
peak whereas the effective scattering rate does not - there 
is no vanishing vertex correction in the effective scatter¬ 
ing rate to control this divergence. This is an indication 
that spin exchange and spin flip processes become very 
important near the critical temperature. 

An interesting consequence of the effective interaction 
used in this paper is that the peak rates of momen¬ 
tum relaxation at T = only indirectly depend on J 
through the temperature dependence in 1/ sinh^(/3hw/2), 
and not through the effective interaction Wq(a;). This 
is most easily understood by approaching the ferromag¬ 
netic transition from the paramagnetic side. The non¬ 
interacting spin response for impurities can be cast as 
Xsjs^(T) = —(T'c/T)[J^S^A^(0)]“^, so that the effective 
interaction reduces to 


fUq(w) 


s2 7v(o)(r/7;) + xiS.(9>^) 


(58) 


At r = Tc, there is neither an explicit dependence on J 
nor on T^, we are left with an interaction that only de¬ 
pends on the non-interacting carrier response. The effect 
of this is l/sinh^(/3Iia;/2) plays the primary role for the 
temperature dependence. In the small frequency limit, 
this factor goes as T^. Then, higher critical tempera¬ 
tures have greater scattering rates associated with them; 
this is verified by Eig. 5a. On the other hand, the heights 
of these scattering peaks are largely independent of the 
carrier density n^, seen in Eig. 5b. The primary contri¬ 
bution of carrier density comes from the cut-off frequency 
hoj < ep'y in each of the response functions. 

Turning our attention to spin relaxation, the expres¬ 
sions derived for DP spin relaxation rely on the assump¬ 
tion that scattering lifetimes are shorter by roughly an 
order of magnitude than spin lifetimes. Else, spins are 
free to process in a spin-orbit field and scattering does not 
appreciably affect the rate of spin relaxation. This is 
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especially relevant for hole spins in III-V semiconductors, 
as they experience stronger spin-orbit interactions than 
electrons."*’^ In Fig. 4, we plot several spin relaxation rates 
for different exchange interaction strengths J. It can be 
quickly verihed by a comparison of timescales that the 
scattering rates in Fig. 5 are indeed fast enough to in¬ 
hibit spin relaxation in the DP mechanism. Near the 
ferromagnetic transition, the mechanism becomes quite 
ineffective. 

An interesting consequence of carrier spins relaxing by 
spin exchange with magnetic impurities is the total spin 
of the system is conserved, which is obviously not the 
case for spin-orbit based spin relaxation. Near the ferro¬ 
magnetic transition, this total-spin-conserving relaxation 
mechanism is very efficient and a quasi-equilibrium spin 
polarization is reached quickly. From here, spin-orbit 
based spin relaxation continues to relax carrier spins at 
a much slower rate until a hnal stage of equilibrium is 
reached. This two-stage relaxation process should be ob¬ 
servable by a kink in the time evolution of the magnetiza¬ 
tion. From simple calculations reported in Appendix B, 
we predict this quasi-equilibrium state to have a magne¬ 
tization in the z-direction of 


^ _ _ yAcrya Asyaj _ 

iXcrOs + Xsgl){X<T + Xs) 




We also studied the relative effectiveness of these spin 
relaxation mechanisms near a ferromagnetic instability. 
In the Dyakonov-Perel mechanism, spin fluctuations act 
as effective scatterers, signihcantly inhibiting the rate of 
spin relaxation by precession in a spin-orbit field. When 
this occurs, spin-flip based spin relaxation mechanisms 
become dominant. When these spin-flips are due to ex¬ 
change with magnetic impurities, the total spin of the 
system is conserved and two stages of spin relaxation are 
expected: fast relaxation occurs due to spin exchange and 
leads to a quasi-equilibrium polarization; this is followed 
by slower relaxation which dissipates spin to a spin-orbit 
field. 
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Appendix A: GKB ansatz for a general spin-indexed 
Green’s function 


where Mz^i is the initial magnetization. The magnetic 
susceptibilities here do not have g-factors absorbed. This 
highlights the importance of the ^-factors in dehning the 
quasi-equilibrium state. 

Further comparing the spin relaxation mechanisms, a 
reasonable question to ask is whether the rates of MI, DP, 
and EY can be summed according to Matthiessen’s rule. 
This would not be appropriate, as the MI mechanism re¬ 
laxes spins to a different equilibrium state than DP and 
EY. Especially considering the timescales at which each 
spin relaxation mechanism operates near the ferromag¬ 
netic transition, it makes more sense to order the events 
as MI spin relaxation to reach a quasi-equilibrium state, 
and then EY -|- DP spin relaxation to bring the system 
to its hnal equilibrium state where the total spin of the 
system has decreased. 

VI. CONCLUSION 

We worked with a Zener model of dilute magnetic semi¬ 
conductors, which exhibits a ferromagnetic transition due 
to itinerant carriers interacting via spin exchange with 
localized magnetic impurities. We used this model to 
derive effective mean-held spin interactions between car¬ 
riers for many-body calculations. From here, we found 
analytic expressions for spin-hip spin relaxation, due to 
exchange of spin between carriers and magnetic impuri¬ 
ties, and Dyakonov-Perel. These expressions are valid for 
degenerate carriers with arbitrary spin polarization and 
can be used with a variety of spin-dependent dynamic 
interactions. 


In the time-domain, the GKB ansatz is written^^ 

= 

(^ 2 ,^ 2 ) — (AI) 

We consider here Green’s functions which are matrices in 
spin-space. The retarded and advanced Green’s functions 
are equilibrium quantities, diagonal in spin-space: 

G'Ll(ti,t2) = T*0[±(ti-t2)]e-*^-(‘-‘=^ (A2) 

When the non-equilibrium time-diagonal Green’s func¬ 
tions are assumed to be diagonal in spin-space, they can 
be considered time-independent under the Markovian ap¬ 
proximation. If off-diagonal elements are expected, there 
remains a time propagator which depends on the energy 
difference between spin states: 

= (A3) 

where p is a density matrix, = p, and p^ = 1 — p. 
Then, the spin-indexed GKB ansatz cannot be simply 
written in terms of ti — t 2 , but instead has the following 
form: 

Gtp{ti,t2) = (A4) 

Appendix B: Quasi-equilibrium magnetization 

We seek the intermediate state of spin relaxation where 
carriers have quickly drained their spin into magnetic im¬ 
purities and begin a slower loss of spin to orbital angular 
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momentum via spin-orbit based spin relaxation mecha¬ 
nisms. Spin-exchange interactions between carriers and 
magnetic impurities preserve the total spin (cr) -|- (S), but 
as these particles in general have different g-factors, the 
magnetization is not preserved. In what follows, we con¬ 
sider the z-components of all quantities without explicitly 
denoting this. 

The equilibrium spin densities are determined by min¬ 
imizing the free energy 

FB{{a),{S)) = 

- hb9s{S)B (B1) 

with respect to (cr) and (S) at constant field B. When 
we remove the magnetic field, we must minimize the free 
energy subject to the constraint (cr) -|- (S) = constant. 
This constant is enforced by a Lagrange multiplier B 
which couples with equal strength to each spin. In other 
words, we have physical coupling —^BB{ga{cr) gs{S)) 
and “Lagrange coupling” —fj.BB{{a) + {S)). The mag¬ 
netic response of the a or S spins to B is exactly the 
same as the magnetic response of these spins to physical 
fields i?cr = B/ga and Bs = B/gs- 

Our system starts slightly out equilibrium so that the 
magnetization is initially 

Mi = gBga{<j)i + gBgs{S)r 

= -{Xa + Xs)B, (B2) 


where we have used gBga{<^)i = —X^B and g.Bgs{S)i = 
—XsB. Spin exchanges between carriers and magnetic 
impurities until equilibrium is reached: 

M^q = gsgai^) eq T 9393 {B) eq^ 

where 9Bga{<j)eq = -X<y{B/ga) and 9Bgs{S)eq = 
~Xs{B / gs). The Lagrange multiplier field B is deter¬ 
mined by the condition 


yielding 


and finally 


a)eq + {S)eq = {<j) i + {S)i 


(B4) 


B = gs9< 


Xa9S + Xsgg 
X<y9l + Xsgl 


B. 


(B5) 


_ ixcrgs +Xs9af „ 

■‘^eg — 9 I 9 

Xags + Xsgi 
^ (x^gs + Xsgaf 

{xc^gl + XsgDixc, + Xs 




(B6) 


^ T. Jungwirth, J. Sinova, J. Masek, J. Kucera, and A. H. 

MacDonald, Rev. Mod. Phys. 78, 809 (2006). 

^ A. H. MacDonald, P. Schiffer, and N. Samarth, Nat. 
Mater. 4, 195 (2005). 

® H. Ohno, H. Munekata, T. Penney, S. von Molnar, and 

L. L. Chang, Phys. Rev. Lett. 68, 2664 (1992). 

^ H. Munekata, H. Ohno, S. von Molnar, A. Segmiiller, L. L. 

Chang, and L. Esaki, Phys. Rev. Lett. 63, 1849 (1989). 

® H. Ohno, A. Shen, F. Matsukura, A. Oiwa, A. Endo, 

S. Katsumoto, and Y. lye, Appl. Phys. Lett. 69, 363 
(1996). 

® H. Ohno, Nat. Mater. 9, 952 (2010). 

T. Diet!, Nat. Mater. 9, 965 (2010). 

® T. Jungwirth, J. Wunderlich, V. Novak, K. Olejm'k, B. L. 
Gallagher, R. P. Campion, K. W. Edmonds, A. W. Rush- 
forth, A. J. Ferguson, and P. Nemec, Rev. Mod. Phys. 86, 
855 (2014). 

® H. V. Lohneysen, A. Rosch, M. Vojta, and P. Wolfle, Rev. 
Mod. Phys. 79, 1015 (2007). 

M. 1. D’yakonov and V. 1. Perel’, Zh. Eksp. Teor. Fiz. 60, 
1954 (1971), [Sov. Phys. JETP 33, 1053 (1971)]. 

Y. Yafet, Phys. Rev. 85, 478 (1952); R. J. Elliott, Phys. 
Rev. 96, 266 (1954). 

M. Wu, J. Jiang, and M. Weng, Phys. Rep. 493, 61 (2010). 
M. Dobrowolska, K. Tivakornsasithorn, X. Liu, J. K. Fur- 
dyna, M. Berciu, K. M. Yu, and W. Walukiewicz, Nat. 
Mater. 5, 444 (2012). 

J. Masek, F. Maca, J. Kudrnovsky, O. Makarovsky, 


L. Eaves, R. P. Campion, K. W. Edmonds, A. W. Rush- 
forth, C. T. Foxon, B. L. Gallagher, V. Novak, J. Sinova, 
and T. Jungwirth, Phys. Rev. Lett. 105, 227202 (2010). 
G. M. Dalpian and S.-H. Wei, Phys. Rev. B 72, 115201 
(2005). 

M. Roever, D.-D. Mai, A. Bedoya-Pinto, J. Malindretos, 
and A. Rizzi, physica status solidi (c) 5, 2352 (2008). 

V. Novak, K. Olejm'k, J. Wunderlich, M. Cukr, 
K. Vyborny, A. W. Rushforth, K. W. Edmonds, R. P. 
Campion, B. L. Gallagher, J. Sinova, and T. Jungwirth, 
Phys. Rev. Lett. 101, 077201 (2008). 

E. H. Hwang and S. Das Sarma, Phys. Rev. B 72, 035210 
(2005). 

F. V. Kyrychenko and G. A. Ullrich, Phys. Rev. B 80, 
205202 (2009). 

F. V. Kyrychenko and G. A. Ullrich, Phys. Rev. B 83, 
205206 (2011). 

T. Dietl, H. Ohno, F. Matsukura, J. Gibert, and D. Fer- 
rand. Science 287, 1019 (2000). 

M. Glazov and E. Ivchenko, JETP 99, 1279 (2004). 

M. Q. Weng and M. W. Wu, Phys. Rev. B 68, 075312 
(2003). 

M. Glazov and E. Ivchenko, JETP Letters 75, 403 (2002). 
M. D. Mower, G. Vignale, and I. V. Tokatly, Phys. Rev. 
B 83, 155205 (2011). 

2® J. H. Jiang, Y. Zhou, T. Korn, G. Schuller, and M. W. 
Wu, Phys. Rev. B 79, 155201 (2009). 

K. Shen and M. W. Wu, Phys. Rev. B 85, 075206 (2012). 






12 


J. H. Jiang and M. W. Wu, Phys. Rev. B 79, 125206 (2009). 
C. Zener, Phys. Rev. 81, 440 (1951). 

C. Kittel, Introduction to Solid State Physics, 8th ed. (John 
Wiley & Sons, Inc., Hoboken, NJ, 2005). 

A. S. Arrott, J. Appl. Phys. 103, 07C715 (2008). 

O. Morandi, P.-A. Hervieux, and G. Manfredi, Phys. Rev. 
B 81, 155309 (2010). 

G. Stefanucci and R. van Leeuwen, Nonequilibrium Many- 
Body Theory of Quantum Systems (Cambridge University 
Press, 2013). 

H. Hang and A.-P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors , Springer Series in Solid- 
State Sciences, Vol. 123 (Springer Berlin Heidelberg, Berlin 
Heidelberg, 2008). 

P. Lipavsky, V. Spicka, and B. Velicky, Phys. Rev. B 34, 
6933 (1986). 

F. Aryasetiawan and O. Gunnarsson, Rep. Prog. Phys. 61, 


237 (1998). 

G. Giuliani and G. Vignale, Quantum theory of the elec¬ 
tron liquid (Cambridge University Press, New York; Cam¬ 
bridge, 2005). 

V. Spicka and P. Lipavsky, Phys. Rev. Lett. 73, 3439 
(1994). 

G. Pikus and A. Titkov, in Optical Orientation, Modern 
Problems in Condensed Matter Sciences, Vol. 8, edited by 

F. Meier and B. Zakharchenya (Elsevier, 1984) pp. 73 - 
131. 

J. Fabian, A. Matos-Abiague, C. Ertler, P. Stano, and 

I. Zutic, Acta physica Slovaca 57, 565 (2007). 

J. Sykes and G. Brooker, Ann. Phys. 56, 1 (1970). 

G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

M. 1. Dyakonov, Spin Physics in Semiconductors 
(Springer, Verlag Berlin Heidelberg, 2008). 

D. J. Hilton and C. L. Tang, Phys. Rev. Lett. 89, 146601 

( 2002 ). 


